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Zinc is an essential mineral for humans and plants and is involved in many physiological 
and biochemical processes. In humans, Zn deficiency has been associated with 
retarded growth and reduction of immune response. In plants, Zn is an essential 
component of more than 300 enzymes including RNA polymerase, alkaline phosphatase, 
alcohol dehydrogenase, Cu/Zn superoxidase dismutase, and carbonic anhydrase. The 
accumulation of Zn in plants involves many genes and characterization of the role of 
these genes will be useful in biofortification. Here we report the identification and 
phlyogenetic and sequence characterization of the 23 members of the ZIP (ZRX IRT like 
protein) family of metal transporters and three transcription factors of the bZIP family in 
Phaseolus vulgaris L. Expression patterns of seven of these genes were characterized in 
two bean genotypes (G19833 and DOR364) under two Zn treatments. Tissue analyzed 
included roots and leaves at vegetative and flowering stages, and pods at 20 days 
after flowering. Four of the genes, PvZIP12, PvZIP13, PvZIPW, and Pv bZIPl, showed 
differential expression based on tissue, Zn treatment, and/or genotype. PvZIP12 and 
PvZIP13 were both more highly expressed in G19833 than DOR364. PvZIP12 was most 
highly expressed in vegetative leaves under the Zn (— ) treatment. PvZIPW was highly 
expressed in leaf tissue, especially leaf tissue at flowering stage grown in the Zn (— ) 
treatment. Pv bZIPl was most highly expressed in leaf and pod tissue. The 23 PvZ\P 
genes and three bZIP genes were mapped on the DOR364 x G19833 linkage map. 
PvZIP12, PvZIP13, and PvZIPW, Pv bZIP2, and Pv bZIP3 were located near QTLs for Zn 
accumulation in the seed. Based on the expression and mapping results, PvZIP12 is a 
good candidate gene for increasing seed Zn concentration and increase understanding 
of the role of ZIP genes in metal uptake, distribution, and accumulation of zinc in 
P. vulgaris. 
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INTRODUCTION 

Dry beans (Phaseolus vulgaris L) are the most highly consumed 
whole food legume in the world. Beans are a food security crop 
for small farmers and urban poor in many African and Latin 
American countries (Siddiq and Uebersax, 2012). In contrast to 
many other staple crops, beans are rich in a variety of nutri- 
ents, including protein, fiber, folate, and minerals (Juliano, 1999). 
Beans are also a good source of dietary iron and zinc. According 
to the USDA Nutrient Database, a 100 g of cooked beans provides 
an average of 2 mg Fe and 1 mg Zn and the Estimated Average 
Requirement for Fe ranges from 3 to 23 mg per day and 2.5- 
10.9 mg per day per Zn depending on age and gender (Ahuja 
et al., 2012). Meeting the Fe and Zn dietary requirements is a 
challenge for many people. An estimated two billion people suffer 
from iron deficiency, which is a major cause of anemia (Rastogi 
and Mathers, 2002; Balarajan et al., 2011). Zinc deficiency is also 
widespread, with an estimated 48% of humans at risk, especially 
populations consuming vegetarian diets rich in unrefined cere- 
als (Sandstead, 1991). In humans, Zn deficiency can be expressed 



through diverse symptoms including reduced immune function, 
fetal brain cell development and child's growth, reproductive, 
and cognitive development (Hambidge, 2000). Biofortification of 
staple foods, including dry beans, with Fe and Zn is one agri- 
cultural based approach being developed and applied to combat 
micronutrient malnutrition (Bouis et al., 2011). While average 
dry bean Fe and Zn levels are 55 mg kg -1 and 34 mg kg -1 
respectively, three-fold genotypic variation in both Fe and Zn lev- 
els exist within the species (Blair et al, 2009 and Islam et al., 
2002). 

This existing variation makes breeding common beans a viable 
biofortification approach. Significant progress has been achieved 
in Fe biofortification of beans through conventional breeding as 
illustrated in the recent release of five high Fe bean varieties in 
Rwanda (Saltzman et al., 2013). Zinc biofortification has lagged 
behind that of Fe-biofortification perhaps because of lower quan- 
tities of Zn in the seeds but also perhaps less incentive because of 
the difficulty in assessing Zn nutritional status in humans. While 
there are biomarkers to asses Fe deficiency readily in humans, 
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no such biomarkers are yet available for Zn, although recently 
a potential biomarker (dematin) has been identified (Ryu et al., 
2012). 

In addition to relying solely on phenotypic selection to 
increase seed Fe and Zn levels, there has been an effort to under- 
stand the genetic control of seed Zn and Fe accumulation. Since 
2009, at least five QTL studies have been published for seed 
micronutrient levels. In total, 38 QTLs were associated with zinc 
accumulation, explaining 15-40% of the variability. These stud- 
ies have been in inter gene pool populations (Blair et al., 2009, 
2010a,b), Andean populations (Cichy et al., 2009; Blair et al., 
2011) and Mesoamerican populations (Blair et al, 2010a). QTL 
studies have yet to be applied to marker assisted selection. There 
has also been limited effort in identifying genes underlying QTL 
for Fe and Zn. Discovery of genes involved in increased seed 
Fe and Zn levels would be useful for biofortification efforts in 
beans and possibly also as targets for transgenic biofortification 
approach in other crops. 

The Zrt and Irt-like protein (ZIP) family is well characterized 
for its role in Zn transport and to a lesser extent it role in Fe trans- 
port (Eide et al., 1996). The ZIP family is well conserved among 
bacteria, fungi, protozoa, animals, and plants (Grotz et al., 1998; 
Chen et al., 2008). ZIP proteins are predicted to have eight trans- 
membrane domains with a histidine motif which may be part of 
an intramembranous heavy metal binding site that plays a role in 
the transport pathway for the minerals that are transferred (Eng 
et al., 1998). ZIP transporters have been implicated in Zn uptake, 
transport of Zn in leaves and translocation to seeds, embryo, 
endosperm, and seed coat (Waters and Sankaran, 201 1). Previous 
information on the role of ZIP genes in Zn movement throughout 
the plant come from expression analysis, yeast complementation 
and Zn hyper accumulator mutants. In A. thaliana 15 members 
have been identified and characterized, revealing a wide variety of 
localization and function (Milner et al., 2012). AfZIPs have been 
detected mainly in the roots, shoots (Milner et al., 2012). In rice, 
17 ZIP coding sequences (CDS) were identified. They have been 
evaluated in roots, shoots, and panicles of efficient and inefficient 
genotypes (Grotz et al., 1998; Guerinot, 2000; Connolly et al., 
2002; Weber et al, 2004; Chen et al., 2008; Shanmugam et al, 
2011; Milner et al., 2012). In Medicago truncatula, six genes were 
identified in roots and leaves which were upregulated under Zn 
deficiency and three of them restored yeast growth on Zn-limited 
media (Lopez-Millan et al., 2004). In Glycine max, GmZIPl has 
been detected in nodules and was highly selective for Zn in a func- 
tional complement in yeast (Moreau et al., 2001). In Vitis vinifera, 
VvZIP3 was expressed in developing flowers and its expression 
was correlated with high Zn accumulation in this tissue (Gainza- 
Cortes et al., 2012 and Afoufa-Bastien et al., 2010). Analysis of this 
family in different species demonstrates the importance of these 
genes in Zn transport. 

Another important gene family related with Zn transport 
is the bZIP family. This family has been well characterized 
in Arabidopsis with 75 members divided in 10 groups based 
on conserved motifs that reflect functional similarities (Jakoby 
et al, 2002). Group F includes bZIP19, bZIP23, and bZIP24. 
These transcription factors contain a DNA binding domain, a 
leucine zipper dimerization motif and histidine -rich motif which 



are essential for responding to low Zn supply in Arabidopsis 
(Assuncao et al., 2003; Assuncao et al., 2010). 

With the recent release of the P. vulgaris genome sequence 
(Phaseolus vulgaris vl.0, DOE-JGI and USDA-NIFA, http://www. 
phytozome.net/commonbean), it is possible to identify candi- 
date genes for seed Fe and Zn levels. Characterization of genes 
related to Zn homeostasis in P. vulgaris will provide useful infor- 
mation on specific target genes in the biofortification breeding 
effort. This research has identified and characterized of 23 mem- 
bers of the PvZIP gene family. Three members of a second family 
of genes, bZIP transcription factors, were also characterized sim- 
ilarly. The relative expression of genes from both the ZIP and 
bZIP families was characterized in various tissues and stages of 
development in two P. vulgaris genotypes, DOR 364 and G19833 
grown under two Zn treatments. Selected ZIP and bZIP genes 
were also located on a linkage map overlaid with QTL locations 
for Zn accumulation in seed. 

MATERIALS AND METHODS 

PLANT MATERIAL AND PHENOTYPIC DATA 

Two bean genotypes were evaluated in this study, DOR364, a 
small seeded, high yielding improved cultivar from the Middle 
American genepool and G19833, a large seeded landrace from the 
Andean genepool known for its tolerance to low P soils (Beebe 
et al., 2006). These genotypes also exhibit contrasting seed min- 
eral levels as shown in field trials in Darien, Colombia. DOR364 
had 49 mg kg" 1 Fe while G19833 had 75.5 mg kg" 1 , and DOR364 
had 21.7 mg kg" 1 Zn while G19833 had 29.9 mg kg" 1 (Blair 
et al, 2009). DOR364 and G19833 were specifically chosen for 
this study because valuable genetic information exists for the 
lines. A recombinant inbred line (RIL) between these parents was 
developed by single seed descent at the International Center for 
Tropical Agriculture (CIAT), Colombia. It consists of 87 individ- 
uals and has a linkage map of 499 single copy markers with a 
coverage of 2306 cM (Galeano et al., 2011). This population has 
been used by different research groups for map saturation and 
QTL identification associated to biotic and abiotic traits (Blair 
et al., 2009) and QTL positions for seed minerals (Blair et al., 
2003, 2009; Beebe et al, 2006; Galeano et al, 201 1). 

IDENTIFICATION OF Pv ZIP AND Pv bZIP GENES AND PHYL0GENETIC 
ANALYSIS 

ZIP genes in P. vulgaris were identified using the sequences of 18 
Arabidopsis thaliana ZIP genes (http://www.arabidopsis.org/). 
The program tBlastn was used to compare the Arabidopsis ZIP 
genes against the bean genome (Phaseolus vulgaris vl.0, DOE-JGI 
and USDA-NIFA, http://www.phytozome.net/commonbean). 
These sequence data were produced by the US Department 
of Energy Joint Genome Institute. Conserved domains in 
each predicted transcript was verified using Pfam 26.0 protein 
database (http://pfam.sanger.ac.uk/) to confirm the reliability 
of the match with the ZIP family. The CDS for each gene was 
aligned with genomic DNA sequence to confirm splice signals 
in boundaries between introns and exons. The P. vulgaris ZIP 
genes were assigned unique names from PvZIPl to PvZIP19 and 
PvIRTl to PvIRT4. These names do not relate to naming of ZIP 
genes in others species. Since this gene family characterization 



Frontiers in Plant Science | Plant Physiology 



July 2013 | Volume 4 | Article 286 | 2 



Astudillo et al. 



ZIP genes in Phaseolus vulgaris 



is based on an incomplete genome sequence, the existence of 
additional ZIP genes in the bean genome is a possibility. 

Three Pv bZIP genes were identified in the P. vulgaris genome 
based on sequences of bZIP19, bZIP23, and BZIP24 reported by 
Assuncao et al. (2010) in Arabidopsis. Identification of the new 
bZIP genes was based on the homology with the Basic Leucine 
Zipper Domain (bZIP domain). 

Sequence alignments, phylogenetic analysis, tree estimation 
using bootstrapping and graphs of each gene were performed 
using ClustalW (Larkin et al., 2007) using the program Geneious® 
6.0.3, created by Biomatters (build 2012-11-06 10:52). 

In silico MAPPING OF Pv ZIP AND Pv bZIP GENES 

Each of the 23 putative ZIP transport protein genes and 3 putative 
bZIP transcription factor genes were mapped in silico to a loca- 
tion on the DOR364 x G19833 linkage map based on sequence 
homology with the P. vulgaris genome. This alignment was con- 
ducted with an MS Excel based program MapSynteny (Fernandez 
et al., 2011). 

GENETIC MAPPING OF SELECT MEMBERS OF THE PvZIP AND Pv bZIP 
FAMILY GENES 

Five ZIP genes were also mapped genetically in the DOR364 x 
G19833 population. These five genes were chosen for genetic 
mapping based on their location near QTLs for seed Fe and Zn 
concentration on chromosomes 1, 3, 6, and 8 (Blair et al., 2009). 
The ZIP genes located in silico in these regions were mapped 
genetically in the full set of RILs of the DOR364 x G19833 pop- 
ulation. These include PvZIP2, PvZIP6, PvZIP8, PvZIPB, and 
PvIRT3. Primers were designed to flank ZIP gene intron sequence 
(Table 1). PCR was conducted on DOR364 and G 19833 as a first 
step to test for polymorphisms. The PCR mix contained 2.0 mM 
Mg, 0.2 [lM dNTPs, and 0.3 [lM of each primer. PCR reactions 
were carried out for 3 min at 95° C, followed by 35 cycles of 30 s 
at 95° C, 30 s at 55 or 60° C (based on the annealing temperature 
of each primer), and a final period of 5 min at 72° C. Products 
were visualized on agarose gels to verify amplification and iden- 
tify insertion/deletions that had potential to serve as molecular 
markers. To increase the possibility of finding polymorphisms for 
monomorphic products, the SSCP technique (from single strand 
conformational polymorphism) was used, which is based on 
detection of conformational differences of single stranded DNA 



fragments due mobility shifts in non-denaturing polyacrylamide 
gel electrophoresis (Orita et al., 1989) such as MDE acrylamide 
gels (MDE Gel Solution 250ML Lonza NJ, USA) as described 
in Galeano et al. (2009). For genetic mapping, Mapdisto soft- 
ware version 1.7 Beta 132 (Lorieux, 2012) was used to locate 
the position of the ZIP genes on the DOR364 x G19833 genetic 
map reported by Galeano et al. (2011). The command place 
locus was used to located the ZIP genes, using as criteria the 
highest LOD value and lowest recombination rate. The posi- 
tion of each ZIP gene was confirmed using the Ripple order 
command. 

QTL DATA AND ANALYSIS 

Phenotypic data for seed Fe and Zn concentration from Popayan 
and Darien Colombia in 1998 and 2003 were reported for this 
population in Blair et al. (2009). Additionally, seed Fe and Zn 
concentration from the same locations in 2006 (not previously 
reported) were used for QTL analysis with the linkage map 
reported in Galeano et al. (201 1). QTL cartographer v. 2.5 (Wang 
et al., 2012) was used to find QTLs following the same parameters 
described in Blair et al. (2009). 

EXPRESSION ANALYSIS OF SELECT PvZIP AND PvbZIP 
Plant growing conditions 

Seeds of DOR364 and G19833 were surface sterilized and planted 
in 500 ml clay pots with 3:1 Sunshine Brand premium grade 
vermiculite (Sunshine Brand, TX, USA) and horticultural grade 
perlite (Industries, Inc., MA, USA). Half strength Hoagland solu- 
tion (3 mM KNO3, 2 mM Ca (N0 3 ) 2 x 4H 2 0, sequestrene DTPA 
10% Fe, l.OmM MgS0 4 x 7H 2 0, 23.1 mM H3B03, 0.38 mM 
ZnS04 x 7H2Q 0.16 mM CuS04 x 5H20, 4.6 mM Mo04 x 
2H20, 1M KH2P04 (pH to 6.0) was applied to pots a rate of 
400 ml three times per week. Two Hoagland solution treatments 
were employed: (1) Zn (+) where Zn was added as ZnSC>4 x 
7H 2 0 and (2) Zn (— ). A total of three pots per genotype were 
planted and each one was designated as a biological replicate. 
The experiment was arranged as a randomized complete block 
design. Plants were grown in a growth chamber (1.86 m 2 ) with 
a photoperiod of 16 h light and 8h dark and an average of tem- 
perature of 29°C/20°C (day/night). For the vegetative samples, 
roots, and leaves were collected from the vegetative 3 stage (V3), 
when the third trifoliate leaf was unfolded at node 5. Leaf and 



Table 1 | Primer list for gene expression analysis via RT-qPCR and genetic mapping. 



Gene Sequence Approach 



PvZIP12 


GGGCAGAGGCAAGTGCAGGG 


GGGCGTGATGGAGATGCAGGA 


RTqPCR 


PvZIP13 


CGCGCTCTTCGATTGCCAGGT 


CCACCGGCGTGTAGTGCGTA 


RTqPCR 


PvZIP13 


GCGGTGGCTCGTTGAGTATT 


TGCTATGAGGTCAACAAGAGCC 


Mapping 


PvZIP16 


TGCACGGTTGATGGCGACGG 


ACGGAACTCCTTCGCCATCGT 


RTqPCR 


PvlRT3 


AGAATAACACCATCCCCAAAATTA 


AGTCACTATGGGAATGTCACAGAA 


RTqPCR 


PvlRT3 


AATGCACATCGTGGGGATGC 


G G CTTTAAACTG C G CTTG G G 


Mapping 


bZIPI 


ATGCAACCCACCTGGCCCTGATGCT 


TGCCTGCCCTTGTAGTTTCCTCGCT 


RTqPCR 


bZIP2 


ATCGGGAGAAGAAGAAGGCTCGCGC 


TCCGGCCCCTTATGTCCACCAGCAA 


RTqPCR 


bZIP3 


GCAGCAGTTCTTGAGCGTGGAGGCT 


TGAAGGTGGTGTTGCCGAAACCTGCA 


RTqPCR 


Pvactinll 


TGCCATCCAGGCCGTTCTTTCA 


GGGGACTGTGTGGCTGACACC 


RTqPCR 
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root samples collected at flowering were harvested at the R2 stage 
when 30% of the flowers were opened. Pod samples were col- 
lected at 20 days after flowering. Plant tissue was collected in 
labeled sterilized tubes of 50 ml in liquid nitrogen and stored 
at -80°C. 

RNA EXTRACTION AND REAL-TIME QUANTITATIVE PCR 

About 2 g of tissue from each sample collected was ground in 
liquid N2. Total RNA from root and leaf tissue of two develop- 
mental stages was extracted by RNeasy Plant Mini Kit (Qiagen). 
Pods were extracted following a protocol optimized for high 
starch samples (Li and Trick, 2005). Total RNA was stored in 
aliquots at — 80° C. The concentration of RNA was quantified 
through Quant-iT™ RiboGreen (Invitrogen). Two u,g of RNA 
of each sample were treated with DNase I and purified by 0.1 
vol of 3M sodium acetate (pH 5.2) and 3 vol of 100% ethanol. 
cDNA synthesis was carried out by High Capacity cDNA Reverse 
Transcription Kits (Applied Biosystem), using 1 |ig of RNA. 
cDNA concentration was measured by Quant-iT™ PicoGreen 
(Invitrogen). 

The relative expression levels of eighth ZIP genes; PvZIP2, 
PvZIP7, PvZIP6, PvZIPU, PvZIP13, PvZIP16, PvZIPW, and 
PvIRT3 and three transcription factors belonging to the bZIP 
family, bZIPl, bZIP2, bZIP3 were measured using RT-qPCR. 
Primers for RT-qPCR were designed for each gene in such a way 
that they spanned one or two exons in genes with intronic regions 
to detect genomic DNA contamination (Table 1). Quantification 
of all transcripts was performed using the Superscript III 
Platinum SYBR Green One-Step qRT-PCR Kit (Invitrogen, 
Carlsbad, CA) according to the manufacturer's instructions. In 
total three technical replicates of 50 ng of cDNA for each biolog- 
ical replicate of all tissues were used as template. Ten-fold serial 
dilutions were used to determine the efficiencies of each primer. 
RT-qPCR master-mix was prepared as follows: 1 [il of diluted 
cDNA, 5 |xl of 2X SYBR Green Reaction Mix, 0.5 u,l 3 pmol of 
each primer and nuclease-free water in a final volume of 10 
The StepOnePlus u Real-Time PCR System (Applied Biosystems) 
was used for amplification and fluorescence measurement of each 
transcript at each temperature step and cycle during the reac- 
tion. Thermal cycling conditions consisted of lOmin at 95° C 
followed by 40 cycles of 1 5 s at 95° C and 45 s at 60° C. The identity 
and purity of the amplified product was checked through analy- 
sis of the melting curve carried out at the end of amplification. 
Relative gene expression was calculated using the comparative CT 
method (Livak and Schmittgen, 2001). bActin was used as a ref- 
erence gene and root in vegetative stage Zn (— ) treatment as a 
calibrator (Wen et al, 2005). Fold changes of greater than 2 was 
used as criteria to determine if genes were differentially expressed. 
Statistical analysis was performed using SAS V 9.3 (SAS Institute 
Inc., NC, USA). A repeated measurement analysis (Proc Mixed) 
was performed. Main effects were tested by AN OVA and a prob- 
ability of P < 0.05 was chosen as the level of significance for the 
statistical test. 

QUANTIFICATION OFZn CONCENTRATIONS IN TISSUE 

Plant tissue from two biological replicates of DOR364 and 
G19833 under two zinc treatments was quantified for Zn 



concentration. Tissue was freeze dried and ground to powder 
using a Geno Grinder 2000 (Spex CertiPrep, Metuchen, NJ) 
and zircon grinding balls. Two grams were sent to A&L Great 
Lakes Labs, Inc. Fort Wayne, IN, for mineral analysis by induced 
coupled plasma spectroscopy. 

RESULTS 

IDENTIFICATION OF ZIP FAMILY MEMBERS AND COMPARISON WITH 
H0M0L0GS IN OTHER SPECIES 

Twenty three sequences, including 19 ZIP and four IRT genes 
were identified in the P. vulgaris genome sequence based on 
similarity to ZIP genes in A. thaliana and/or Medicago trun- 
catula. All new genes have full-length CDS containing open 
reading frames (ORF) ranging from 153 to 655 amino acids 
in length. Sequences identified were confirmed in the PFAM 
database based on ZIP transmembrane domain and had E-values 
higher than —10. Peptide sequences of all new ZIP genes iden- 
tified in common bean were aligned with 18 ZIP genes reported 
in the A. thaliana, and M. truncatula (Table 2). A phylogenetic 
neighbor joining tree shows the relationship among ZIP genes in 
P. vulgaris, A. thaliana, andM. truncatula (Figure 1). Alignments 
at the amino acid level predicted eight highly conserved trans- 
membrane domains (Figure 2) and a potential metal binding 
motif containing histidine residues implicated in metal trans- 
port which are highly conserved throughout the entire family 
(Guerinot, 2000; Lopez-Millan et al., 2004). All ZIP genes con- 
tained a histidine motif between transmembrane domain III and 
IV except PvZIP6, PvZIP7, and PvZIP18. The ZIP gene fam- 
ily members in P. vulgaris shared 3-81.4% homology to each 
other. Of all ZIP genes found PvIRT3 was the most closely 
related to Arabidopsis, sharing 59 and 57.3% similarity with genes 
AtIRT3_ATlG60960.1 and AtZIP4_ATlG10970, respectively. 
PvZIPU also showed high similarity with AtZIP6_AT3G30080.1 
at 53.8%. 

Gene structure analysis of ZIP genes in P. vulgaris revealed 
that the 23 genes have different intron-exon structures with a 
wide range of lengths. PvZIP2, PvZIP6, PvZIP7, PvZIP15, and 
PvIRTl, are composed of three exons and two introns. PvZIP3, 
PvZIP5, PvZIP9, PvZIPIO, PvZIPU, PvZIPU, PvZIP16, PvZIP19, 
and PvIRT2 each have four exons and three introns. PvZIP17 and 
PvIRT3 have five exons and four introns. Seven exons were iden- 
tified in PvZIPl and PvIRT4. Many exons (10-14) were present in 
PvZIP4, PvZIP8, PvZIPU, PvZIPU, and PvZIPW. 

Given the importance of some members of the bZIP gene 
family in the regulation of ZIP genes and in turn plant Zn home- 
ostasis, their sequences were also characterized in the P. vulgaris 
genome. The common bean genes bZIPl, bZIP2 and bZIP3 were 
261, 266 and 154 amino acids long, respectively. None of the bZIP 
genes contained introns. The three amino acids sequences encod- 
ing the bZIP genes shared 4.0-38.5% similarity among each other 
and 15-55.4% similarity with bZIP19, bZIP23, and bZIP24 genes 
described in A. thaliana (Assuncao et al., 2010). 

MAPPING OF PvZIP GENES AND QTL FOR SEED Fe AND Zn 
CONCENTRATION 

ZIP and bZIP were mapped in silico on the DOR364 x G19833 
genetic map by aligning ZIP gene sequences and molecular 



Frontiers in Plant Science [ Plant Physiology 



July 2013 | Volume 4 | Article 286 | 4 



Astudillo et al. 



ZIP genes in Phaseolus vulgaris 



o o r- 

CD CD CO LO ^ 

CD T- T- ^ T- 



CM LO CM CM 



r-- co 



OCMr^OCOCOCOCDCMCO 



o o 

LO LO 

o o 

CM CM 

CO CO 

o o 

CD CD 

CO CO 

T3 T3 



LO t3 t3 CD 
O "O "D LD 
CM CD CD CO 



LO LO LO i 

CD CD CD CD 



._ co co co co 

N C C C C 

0 ro ro co ca 

_^ £3 £3 £3 £ 

TT tj tj o o 



INI M M N N 



O O 

LO LO 

CM CM "v 5 

00 CD -h 

O O "O 



O <D CD 



CM O O O 
00 CO CD CO 
O LO CD O 



CO CD CD 



O O O 

CD CD Q 
CM LO " 



CD 0 0 CD 



I*-. O 5f" CD 



\J vl 'k_ CD CO y 
■<-■ LO H D) 

"o "o "d >— 4=; r-~ 



LO LO LO 



t t 



o o ^ ^ 

ro H to CO CO 

LO O CO Uj Uj 

O 3 § o o 

O) O D) ^ 2 

e ™ 2 ^ € 



? 1 1 

I CD CD 



CD CD 



Q.Q.Q.Q.Q.Q.Q.Q.O Q.Q.Q. 

cococococotototoQ-cototo 

CD CD CD CD CD CD CD 



U CJ CJ CJ 



H ± h H H ro fa 
O O CJ * ' 



o o o 



MNNNMMNINININMNINIfNlMN 



Q. ■— O 



CN O CD S fS ^ 
CO CO CM i LO r— 



COO^l-CDCOCMCMCO 



S CM CO CM 
C^J O CM i— 



LLJ <>i ^ 
LO CM CD 



r-- cm 



p >. -P 



I I t e E E 

cu CD o o ro CD 
t t Q. Q. 4— 



+zl +zi CD 
O O CD £ 



•S" "St 
Q_ Q_ 

|SI |SI 



Q. ■- 

CU CD 

J= E 

B 9= 

N < 



S CD CD 



N. M. INI, 



Q_ Q_ 
N N 



CD CD CD 



E 

Q_ 



CO CO CD CD O O 



O Q. "5 c m 



co co co 



CD i3 CD CO 



£5 E 



INI. INI. INI. 



CO CD CD N CD 

9= E 9= J E 



CD 
<Z 

.9 o 

_ Q. 

CD CO 

I § 

Q_ O 

8 o 1 



ro ro 
I — h- 
tr dc 



CO CD "vf 
CM i- -st 
Q_ Q_ Q_ 



M M Nl 



CD ^ CM i— 

^ O CD" to' °. 

CM i^- CO CO LO 

5 h CO N 

^ CD CD CO" & 





cn 


CO 




m 


CO 


CD 






03 


03 


cn 


co" 




LD~ 


LO 

CO 


CM 


CO 










o 






I-." 


o 


CO 


co 


co 


CN 



CD CD CO 

co" o" ^ 



CD 


C50 
O 


05 

cn 

CD 


■* 


CN 


LO 


cn 


cn 






LO 


CO 


CO 


LO 
CO 


,817, 


CD 




o> 




LO 


LO 


cn 



CO CO CO LO 
O O) • 



CM CO CD 



co ^f- cd 



'-r--cooococo'-coo 



- CD CD CD CM CM 



r cm ^ 









CN 


LO 


LO 


LO 


LO 


CD 


CD 


CD 


CD 


CD 


CD 


03 


03 


o 




CN 


oo 


oo 


cn 


cn 


LO 






o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


O 


O 


O 


O 


O 


O 






O 


o 


o 


o 


o 


o 












-c: 


-c: 


-c: 


-c: 




-c: 


-c: 


-c: 


JZ 


JZ 


-C 


-C 


-C 




JZ 


-C 






JZ 


JZ 


JZ 


JZ 


JZ 


O 


O 


O 


O 


O 


O 


O 


O 


O 


U 


U 


O 


O 


O 


O 


O 


O 


O 


O 


o 


o 


O 


O 


O 


O 


O 





















O 




CN 


CO 




LO 


CD 




03 


cn 


















CN 


CO 




LO 


CD 




03 


OT 
























CN 


CO 












Q_ 


Q_ 




Q_ 


Q_ 


Q_ 


Q_ 


Q_ 


Q_ 


Q_ 




CL 


CL 


CL 


CL 


CL 


CL 


CL 


h— 


I — 


I — 


I — 




CN 


00 


ISI 


N 


N 


N 


N 


N 


N 


N 


N 


N 


N 


N 


N 


M 


N 


N 


N 


N 


M 


CC 


CC 


CC 


CC 


CL 


CL 


CL 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 




hJ 


M 


Q_ 


Q_ 


Q_ 


Q_ 


Q_ 


Q_ 


Q_ 


CL- 


Q_ 


Q_ 


Q_ 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


CL 


-Q 


-Q 





E E 

CN LO 

co 

o LO 

5 5 



E E 



b co 

CN OT 



O O O O 



E E 

^ CN 

i- O 
CD 

LO OT 



E E 

OT 

03 CD 

CD CD 

OT OT 

CN CN 

O O 



E EE E 



CN O O 



E E 

LO CO 

OT <- 

O *~~ 

00 CN 



E E 



E E 
2 lo 



O O CN O 



E E 

LO CO 



E E E E E E E 

CN CD CD LO go O 5 

00 'J 

LO LO 



CO £ >n CO 03 CO 

" ™ w a a co 



O O ^ ^ i- 



oooooooooooooooooo 



CLCLCLCLCLCLCLCLCLCLCLQ_Q_CLCLCLCLQ_CLCLCLCLCLCLQ_D_ 



www.frontiersin.org 



July 2013 | Volume 4 | Article 286 | 5 



Astudillo et al. 



ZIP genes in Phaseolus vulgaris 






MtZIP3 

PvZIP2 

PvZIPIO 

PvZIP1 1 

PvZIP9 

PvZIP16 

PvZIP5 

PvlRT3 

PvZIPI 

WWZIP14 

PvZIPI 3 

PvZIPI 5 

AtZIP6 

MtZIP11 

MtZIP5 

MtZIP18 

AtlRT3 

AtZIP4 

AtZIP9 

MtZIPIO 

PvZIP6 

PvZIP7 

MtZIP4 

PvZIP4 

AtJARI 

PvZIPI 2 

AtZIP1 1 

AtZIP2 

MtZIP7 

MtZIP8 

PvZIP8 

PvZIPI 8 

AtZIP29 

MtZIP9 

PvZIPI 4 

PvZIPI 7 

AT3G08650.2 

MtZIP19 

MtZIP17 

WWZIP16 

MtZIP6 

PvlRTI 

PvlRT2 

PvlRT4 

PvZIPI 9 

AtZIPIO 

AtlRTI 

AtlRT2 

AtZIP8 

AtZIP7 

MEIP15 

AtZIP12 

MtZIP12 

MtZIP13 

PvZIP3 

MtZIPI 

AtZIPI 

AtZIP5 

AtZIP3 



0.9 

FIGURE 1 | Phylogenetic tree of homologs ZRT, IRT-like protein family 
in Phaseolus vulgaris, Arabidopsis thaliana, and Medicago truncatula. 

Analysis was based on alignment of amino acid sequences using 
Geneious program v. 6.0.3 and N-J trees were generated. Arabidopsis 



genes are indicated with the ZIP and IRT number used on TAIR 
database. ZIP names used in Medicago were according to Lopez-Millan 
et al. (2004) (ZIP1 to ZIP7). ZIP8 in front were assigned with a 
consecutive number. 



marker sequence in DOR 364 x G19833 against the P. vulgaris 
genome sequence. The results of the in silico mapping indicate 
that ZIP genes are distributed on all P. vulgaris chromosomes 
except 4 and 7. There was a tendency for ZIP genes to cluster 
together, most notably on chromosome 5 and 6 (Figure 3). 

Through in silico mapping, bZIPl was located on chro- 
mosome 5 and bZIP2 and bZIP3 were located near each 
other on chromosome 11. The bZIP3 gene location was 
expected based on the position and sequence of the SNP 
marker g785 which contains the bZIP domain. This marker 



was described as PvMcCleanNDSU2007_ll_g785 (http://cmap. 
comparative-legumes.org). Selected ZIP genes were also mapped 
genetically via DNA polymorphisms in the DOR364 x G19833 
population using the P. vulgaris reference genetic map published 
by Galeano et al. (2011). This map consists of 499 single copy 
markers and 2306 cM of coverage (Figure 3). PvZIP2, PvZIP6 — 
PvZIP8, PvZIP13, and PvIRT3 mapped to chromosomes 1, 5, 6, 
and 9, respectively. 

Once gene markers were mapped, QTLs for seed Fe and Zn 
were also identified on this map. These QTLs include previously 
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sequence located in the variable region between transmembrane domains III 
and IV and fully conserved histidine motifs are indicated by gray lines. The 
eight domains are shown as a red line above the sequences. 



published data for two sites (Blair et al., 2009) as well as QTL 
identified in whole and cotyledon seed mineral evaluation from 
a 2006 planting of the same population in Darien, Colombia. 
QTL analysis in the 2006 evaluation identified new QTLs for Zn 
concentration on chromosomes 1 and 2 and also confirmed the 
QTLs identified by Blair et al. (2009) (Table 3). For seed Fe con- 
centration, 13 QTLs were found on chromosomes 2, 3, 6, 8, and 
11. For seed Zn concentration 11 QTLs were found on 1, 2, 3, 
6, 9, and 11 (Table 3). ZIP genes mapped on chromosomes 3, 
6, 8, and 11 mapped within the regions of QTL for Fe and/or 
Zn. On chromosome 11 two bZIP genes (with genomic posi- 
tion 3,134,797 and 3,709,270 bp) were mapped in silico within 
the region of two QTLs for seed Fe and one QTL for seed Zn 
(Figure 3). Two iMRT genes are present on chromosome 3 (at 
49,001,506 and 49,013,793 bp) and three QTLs for seed Fe con- 
centration mapped between the QTLs (Figure 3). Table 3 shows 
specifically which ZIP genes are located within or nearby QTL for 
seed mineral concentration. 

EXPRESSION ANALYSIS OF PvZIP GENES 

Studies in Arabidopsis, Glycine, Vitis and Medicago indicate that 
ZIP genes may be expressed in roots, leaves, and reproductive 
tissue (Grotz et al, 1998; Lopez-Millan et al, 2004). Many stud- 
ies so far have focused on expression in roots and shoots (Grotz 
et al., 1998; Lopez-Millan et al, 2004; Milner et al., 2012). From 
the perspective of biofortification, it is necessary for a bean plant 
not only to efficiently take up Zn from the soil, but also trans- 
port and accumulate it in vegetative tissue, pods, and ultimately 
seeds. In order to determine the expression profile of members of 
ZIP family and their relevance during the development of com- 
mon bean, relative expression levels were measured by RT qPCR. 
PvZIP2, PvZIP7, PvZIP6, PvZIP12, PvZIP13, PvZIP16, PvZIP18, 
and PvIRT3 genes were selected for this analysis based on their 
location in the genome in relation to presence of QTLs for Zn 
and Fe in the DOR364 x G19833 population. Three tissue types 



were analyzed for gene expression in DOR364 and G19833: roots, 
leaves, and pods. Roots and leaves were collected at two time 
points, one during vegetative growth, and one during flower- 
ing. Pods were sampled 20 days after flowering. Each tissue type 
was selected from plants with two Zn treatments. At four weeks 
after planting, DOR364 and G19833 plants in the Zn (-) treat- 
ment exhibited some Zn deficiency symptoms such as interveinal 
chlorosis, bronzing, and shortening of the internode (Brown and 
Leggett, 1967). The ZIP genes PvZIP12, PvZIP13, PvZIP16, and 
PvIRT3 were expressed in all tissue analyzed (Figure 4). However, 
PvZIP2, PvZIP6, PvZIP7, and PvZIP 18 were undetectable under 
RT qPCR in all tissue types. This finding is also supported by 
pod transcriptome data which also found low to no expression 
for PvZIP2, PvZIP6, PvZIP7, and PvZIP18 (Astudillo et al, in 
preparation). 

Gene expression of PvZIP 12 and PvZIP 16 was induced upon 
Zn status in leaf tissue. PvZIP12 was most highly expressed in 
leaves at vegetative stage under Zn (— ) treatment, especially 
in G19833 (Figure 4). For PvZIP13, G19833 was more highly 
expressed than DOR364 especially in leaves sampled during flow- 
ering. Of each of the ZIP genes studied, PvZIP 16 showed the 
highest differential expression based on tissue type, stages, and 
zinc treatment but not between genotypes. It was 139-848-fold 
more expressed in the leaves than the roots for both genotypes 
and developmental stages. 

EXPRESSION ANALYSIS OF THREE TRANSCRIPTION FACTORS bZIP 

RNA from the same samples described above were also used to 
determine the relative expression of three transcription factors 
Pv bZIPl, Pv bZIP2 and Pv bZIP3, which are homologous to 
Arabidopsis bZIP genes in the zinc homeostasis network (Table 2). 
The common bean homologue bZIPl was detected in roots, 
leaves, (at vegetative stages) and pods but expression pattern did 
not change based on stages and Zn treatment. This gene was more 
highly expressed in leaf tissue sampled during flowering than 
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FIGURE 4 | Relative expression level of PvZIP gene transporters 
and three bZIP transcription factors in genotypes Dor364 and 
G19833 in different tissues and two Zn treatment, (i) roots at 
vegetative stage (V_ROOT— and V_ROOT+), (ii) roots at flowering 



stage (F_ROOT— and F_ROOT+); (iii) leaves at vegetative stage 
(V_LEAF— AND V_LEAF+) stage; (iv) leaves at flowering stage 
(F_LEAF- and F_LEAF+); and (v) pods (POD- and POD+) of plants 
under Zn (— ) and Zn (+) treatment. 
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vegetative tissue in both G19833 and DOR364. Transcripts of Pv 
bZIP2 were detected in roots, leaves, and pods and its expression 
pattern was not influenced by tissue type, developmental stage or 
Zn treatment. Pv bZIP3 was expressed in roots during vegetative 
and flowering with a barely expression in leaves and pods. 

TISSUE ZINC CONCENTRATION 

Zn concentration was determined for DOR364 and G19833 for 
all tissues, developmental stages, and Zn treatments (Figure 5). 
The highest Zn concentration was observed in roots for both 
genotypes and no significant differences were observed in leaves. 
Significant tissue x Zn treatment x genotype interaction was 
found. Plants grown under the Zn (+) treatment had higher levels 
of Zn in pods and seeds than those grown under the Zn (— ) treat- 
ment. Seed Zn levels were 26 and 53% less in the Zn (— ) treatment 
in DOR364 and G19833, respectively. G19833 had higher seed Zn 
levels that DOR364 under the Zn (+) treatment but not under 
the Zn (— ) treatment (Figure 5). 

DISCUSSION 

Common bean is becoming an alternative to dietary supplements 
as a way to improve human health. ZIP metal transporters are 
one of the most important gene families for Zn and Fe cellular 
uptake and translocation in plants (Guerinot, 2000; Chen et al., 
2008; Wu et al, 2009; Adams et al., 2012). Identification of ZIP 
members in P. vulgaris and characterization of their expression 
patterns is useful to increase the understanding of uptake, trans- 
portation, and storage of Zn. This study is a unique combination 
of gene family characterization with physical and genetic map- 
ping and functional expression data that has utility in common 
bean improvement. 

Twenty-three ZIP genes were identified in the P. vulgaris 
genome and genes were annotated and characterized based on 
similarity to other ZIP family members in A. thaliana and M. 



truncatula. According to the total number ZIP family members 
across species the family origin may be from a common ancestor 
that has undergone sequence duplication followed by divergence 
events (D'Ovidio et al., 2004). PvZIP genes clustered on chro- 
mosomes 3, 5, 6, and 9 showed high sequence similarity. The 
close proximity and sequence similarity of many of the ZIP gene 
family members might suggest of gene duplication followed by 
diversification (Yang et al., 2009). On the other hand, hetero- 
geneity in structure and expression in each PvZIP gene suggests 
diversity in functionality. Four of the eight genes evaluated were 
not expressed in none of the tissue analyzed. Loss of function 
in these proteins might be compensated for by duplicate genes 
(D'Ovidio et al., 2004). 

It is important to consider the link between functional varia- 
tion and gene structural differences among ZIP family members. 
In many cases, Zn interacts with cysteines and histidines in pro- 
teins and may determine the ionic selectivity of ion transporters 
(Ramesh et al., 2003; Lopez-Millan et al, 2004). The motif of his- 
tidine in variable regions between transmembrane domain III and 
IV in many ZIPs has been postulated to serve as a potential metal 
ion binding site (Eide et al., 1996; Zhao and Eide, 1996; Grotz 
et al., 1998). For PvZIPs identified in this study, all contained this 
motif except PvZIP6, PvZIP7, and PvZIP 18, interestingly these 
ZIP genes were also not expressed in all tissue analyzed, suggesting 
without the motif they are not functional. 

In Arabidopsis, ZIP genes have been shown to regulate and 
also contribute to the uptake, transport and accumulation of Zn 
(Grotz et al., 1998; Weber et al, 2004; Talke et al, 2006; Lin et al, 
2009; Milner et al., 2013). Here we used RT-qPCR approach to 
obtain a picture of ZIP gene transcription in roots and leaves at 
vegetative and flowering stages, and pods at 20 days after flow- 
ering in P. vulgaris. Some processes such as Zn uptake, have 
been studied in detail, while others such as remobilization of Zn 
from vegetative to reproductive tissues are less well understood 
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FIGURE 5 | Zinc concentration in DOR364 and G19833. Zn concentration V_LEAF+) stage; (iv) leaves at flowering stage (F_LEAF— and F_LEAF+); (v) 
(ppm) in (i) roots at vegetative stage (V_ROOT- and V_ROOT+), (ii) roots at pods (POD- and POD+) and seeds (SEED- and SEED+) of plants under Zn 
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(Gene et al., 2006). The evaluation of gene expression patterns 
based on tissue, Zn treatments, and genotype not only provides 
information on the functionality of the ZIP family genes but also 
may help explain genotypic differences in seed Zn accumula- 
tion. These data indicate differential gene regulation associated 
with the nutritional requirements and possible mechanism of 
partitioning of Zn along the plant. According to ZIP genes char- 
acterization in Arabidopsis approximately half of the genes char- 
acterized are induced in response to Zn deficiency (Grotz et al., 
1998; Talke et al., 2006). 

ZIP gene expression differences in P. vulgaris were related 
to Zn treatments, genotype, and tissue type. Genotypic differ- 
ences in Zn translocation capacity in different organs may be 
an important factor in Zn accumulation in seeds (Hacisalihoglu 
et al., 2004). Observed differences between genotypes could also 
be due to genetic differences and diversity among Andean and 
Mesoamerican gene pools (Blair et al., 2009). Similarly to previ- 
ous studies, G19833 had higher seed Zn level than DOR364 (Blair 
et al., 2009). However DOR364 had higher Zn in its roots as com- 
pared to G19833 suggesting that G19833 can translocate more Zn 
from roots to seeds. 

Zinc plays a specific role in fertilization and pollen grains 
contain very high concentrations of Zn (Fageria et al., 2011). 
At flowering most of the Zn taken up is incorporated into the 
developed seed (liang et al., 2008) so genes highly expressed at 
flowering and in pods such as PvZIP12, PvZIP16, and bZIPl could 
be directly related to Zn remobilization to seeds. 

Although leaves are known as the major source of remobilized 
micronutrients in common bean (Sekara et al., 2005) in rice stems 
are the major source of Zn in the seed (Waters and Sankaran, 
2011). With this study it was not possible to determine how much 
and the source of Zn remobilization. Future studies with radio 
labeled Zn would be warranted to asses Zn remobilization. 

Based on the relative expression values established via RT- 
qPCR, the high Zn concentration in roots did not reflect 
expression values for the ZIP genes evaluated in this tissue. In 
Arabidopsis at least 10 different members of the ZIP family play 
a role in zinc uptake in roots, including ZIP1, 2, 3, 4, 5, 9, 10, 
11, 12, and IRT3 (Van De Mortel et al., 2006). We evaluated four 
of their respective homologous in P. vulgaris and found that they 
were only weakly expressed in roots. 

The DOR364 x G19833 RILs mapping population consists 
of 86 individuals, which are adequate for identifying QTL with 
moderately large effects based on QTLs previously detected (Blair 
et al., 2009, 2011; Galeano et al, 2011). In silico mapping of 
ZIP genes was a successful strategy to locate PvZIP genes aligned 



with QTL for seed Fe and Zn in the bean genome. QTL analy- 
sis was carried out in the current reference bean map (Galeano 
et al., 2011). It is worth noting that on chromosomes 2 and 6 
where PvZIP4 and PvZIP12 and PvZIP13 are located, there are 
QTL for seed Zn concentration. For Fe, the IRT genes are con- 
sidered to be the main transporters for high-affinity iron uptake 
in roots in Arabidopsis (Connolly et al., 2002; Henriques et al., 
2002; Lin et al, 2009). In this study, PvIRTl and PvIRT2 were 
located on chromosome 3 within an important QTL region asso- 
ciated with seed Fe concentration. The Pv bZIP2 and Pv bZIP3 
genes were located on chromosome 1 1 and aligned with the most 
important QTL for Fe and Zn reported in P. vulgaris. There are 
no obvious genotypic differences in expression of these genes 
in G19833 and DOR 364, however. The QTL in this region 
has been found in at least three mapping populations, includ- 
ing Mesoamerican and Andean intra and inter genepool crosses 
(Blair et al., 2009, 2010a,b, 2011). The bZIP transcription fac- 
tors analyzed correspond to genes in Arabidopsis responsible for 
response and adaptation to low Zn supply. In general, PvZIP, 
iMRT, and Pv bZIPs co-localization with QTLs for Fe and Zn 
levels suggest that their function is important in Fe and Zn 
homeostasis in P. vulgaris. In Arabidopsis, the bZIP transcrip- 
tion factors that interacted with ZIP genes were found directly 
upstream of the ZIP genes (Assuncao et al., 2010). In the case 
of P. vulgaris none of the bZIP genes were adjacent to ZIP 
genes. 

This study is the first to characterize the ZIP gene family, report 
the expression profile in various tissues with two genotypes and 
fertilization treatments. It provides evidence of the relationship 
among level of transcripts and QTLs in dry bean seed as was 
identified in PvZIP12 gene. This contribution will be particu- 
larly useful for advancing bean breeding programs. The use of 
such gene marker encoding proteins associated with transport 
and accumulation of Zn and Fe could increase the efficiency and 
accuracy in the selection of breeding materials for biofortification. 
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